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高 超声 速 风 洞 马赫 数 4. 5 喷 管 气动 
设计 和 数值 验证 


黄 风 , 杨 永 能 , 悄 继 斌 , 杨 海 滨 , 张 伟 
(中 国 空气 动力 研究 与 发 展 中 心 ,621000 绵阳 ) 


摘 要 :常规 跨 超声 速 风 洞 进行 马赫 数 4.5 试验 时 ,常常 伴 有 空气 液化 现象 ,造成 试验 数据 可 信 度 
低 , 在 高 超声 速 风 洞 研制 马赫 数 4.5 喷 管 ,具有 对 气流 加 热 的 能 力 ,可 以 提供 更 加 准确 的 试验 数据 。 
目前 国内 0.5 m 量 级 高 超声 速 风 洞 还 不 具备 马赫 数 4.5 的 试验 能 力 。 通 过 无 秋 流 计算 方法 计算 轴 
对 称 喷 管 型 面 ,并 采用 Sivells-Payne 方法 进行 附 面 层 修正 ,然后 进行 数值 验证 ,证 明了 计算 出 的 型 面 
满足 国 军 标 对 马赫 数 的 设计 要 求 , 可 以 投入 加 工 生产 。 
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Aerodynamic design and numerical validation of Mach 4.5S nozzle 
in hypersonic wind tunnel 


HUANG Ju,YANG Yongneng, XU Jibin, YANG Haibin,ZHANG Wei 


(China Aerodynamics Research and Development Center,621000 Mianyang, China) 


jAbstract: The Mach 4.5 test in a conventional tran-supersonic wind tunnel is often accompanied by the 
phenomenon of air liquefaction , resulting in low reliability of the test data. The Mach 4.5 nozzle developed 
in a hypersonic wind tunnel has the ability to heat the airflow ,which can provide more accurate test data. 
At present ,the test capability of Mach 4.5 is not available in China for the 0.5-meter hypersonic wind 
tunnel. The axisymmetric nozzle profile was calculated by the inviscid flow calculation method, and the 
boundary layer was modified by the Sivells-Payne method. Then ,numerical simulation was carried out. It 
was proved that the calculated profile met the GJB design requirements of Mach number and can be put in- 
to production. 
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超声 速 飞行 咒 是 我 国正 在 研制 的 重要 武器 装 。 再 入 段 的 关键 马赫 数 , 也 是 亚 燃 冲压 向 超 燃 冲压 过 
备 ,由 冲压 发 动机 了 驱动 高超 声速 巡航 ,具有 飞行 速 ” 渡 马 赫 数 ,马赫 数 4.5 的 试验 数据 对 飞行 器 研制 至 
度 快 .机 动 飞行 难以 拦截 等 优势 。 马 赫 数 4.5 是 高 。” 关 重 要 。 


超声 速 飞 行 器 .战略 战术 导弹 等 飞行 器 在 上 升 段 和 常规 风 洞 一 般 把 马赫 数 4.5 喷 管 设计 在 超声 速 
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风 洞 ,由 于 超声 速 风 洞 不 具备 气流 加 热 条 件 , 而 马赫 
数 4.5 试验 时 空气 液化 ,造成 试验 数据 可 信和 度 低 。 
目前 国内 工业 性 质 的 0.5 m 量 级 高 超声 速 风 洞 均 无 
马赫 数 4.5 的 风 洞 试验 能 力 。 在 高 超声 速 风 洞 增加 
马赫 数 4.5 喷 管 ,建成 后 可 以 形成 马赫 数 4.5 风 洞 
试验 能 力 , 既 能 满足 当前 武器 装备 型 号 研制 对 马赫 
数 4.5 的 风 洞 试验 能 力 需求 ,也 可 解决 亚 跨 超声 速 
风 洞 与 高 超声 速 风 洞 之 间 试 验 数据 衔接 问题 ,为 其 
提供 试验 设备 和 试验 手段 ,满足 型 号 研制 部 门 从 跨 
超 到 高 超 风 洞 试验 数据 的 需求 。 

本 研究 通过 对 马赫 数 4. 5 喷 管 的 气动 设计 ,并 
对 计算 的 内 型 面 进行 了 数值 验证 ,得 出 了 满足 国 军 
标 要 求 的 马赫 数 喷 管内 型 面 数据 ,可 以 用 于 加 工 生 
i 实物 。 


3 内 型 面 设计 


到 轴 对 称 喷 管 气动 设计 计算 方法 


思 旬 对 称 喷 管 气动 设计 主要 包括 2 部 分 :位 流 型 
天 pi 过 机 和 让 省 音速 段 ) 和 超 音 速 位 流 型 面 边界 
属 哆 正 。 位 流 又 分 为 3 个 区 域 :第 一 区 域 是 亚 音速 
区 第 二 个 区 域 是 由 哈 道 到 拐点 的 喉 道 区 ;第 三 个 区 


change area G-A 


Current source (FGAB) 
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域 是 由 拐点 到 喷 管 出 口 的 转换 区 ( 消 波 区 )。 

在 喉 道 ,外 形 由 单调 递增 的 三 次 曲线 给 出 。 在 
转换 区 ,外 形 主要 由 Cresci 方法 计算 ,如 图 1 所 示 ， 
假定 流动 从 0" 源 点 出 发 , 沿 径 向 流 一 直 膨 胀 到 波 前 
4B, 在 拐点 4 处 具有 最 大 膨胀 角 一 半 锥 角 0 的 锥 形 
流 。CE 下 游 是 均匀 平行 流 区 ,在 ABCE 区 域 中 流动 
由 径 向 流 过 渡 到 平行 流 ,利用 波 前 4B8 和 BC 上 的 已 
知 条 件 ,用 轴 对 称 特征 线 网 格 计算 。 最 后 得 到 7E 位 
流 型 面 。 

在 给 定 0 的 值 ,如 果 B、C 重合 ,这 样 得 到 最 短 
喷 管 。 不 过 在 这 样 的 喷 管 中 , 沿 轴 向 马赫 数 梯度 在 
B 点 不 连续 ,其 曲率 半径 在 拐点 处 不 连续 ,理论 上 ， 
曲率 不 连续 会 导致 在 做 边界 层 修正 时 喷 管 的 边界 层 
为 负 ; 如 曲率 是 连续 的 ,边界 层 的 增长 也 就 是 光 

滑 的 。 

况且 ,即使 边界 层 增长 被 准确 计算 出 来 ,其 喷 管 
型 面 的 黏 性 修正 也 只 是 针对 设计 运转 条 件 的 ,那么 
修正 量 也 只 有 在 计 运 转 条 件 下 才能 够 应 用 到 非 黏 性 


外 形 上 去 ,所 以 离开 设计 运转 条 件 ,边界 层 增 长 的 偏 
差 将 产生 不 连续 ,这些 不 连续 将 顺 着 喷 管 因 机 械 加 
工 不 光滑 的 地 方向 下 移动 ,从 而 导致 试验 段 的 流动 
中 存在 扰动 。 奉 外 形 有 连续 的 曲率 ,3 
条 件 就 可 以 减弱 一 些 。 


这 些 设计 运转 


Uniform 


图 1 轴 对 称 喷 管 气动 设计 计算 图 


Fig.1 Calculation diagram of axisymmetric nozzle aerodynamic design 


为 了 克服 以 上 缺点 ,引入 部 分 消 波 区 4BCD , 喷 
管 增长 了 ,曲率 连续 可 以 满足 ， a i 
在 B 点 马赫 数 Ms 是 可 以 任意 给 定 的 (但 要 比 C 上 
低 ) ,BC 之 间 的 导数 分 布 由 三 。 出 ,方程 中 
系数 由 连续 条 件 确定 。 

喉 道 区 外 形 用 三 次 曲线 计算 , 它 是 单调 递增 到 


具有 0, 的 切 点 上 ,在 4 点 上 二 阶 导数 为 零 ,在 喉 道 处 
满足 连续 条 件 。 


亚 音 速 收缩 段 是 将 稳定 段 来 的 气流 均匀 加 速 至 
声速 。 根 据 高 超声 速 喷 管 的 设计 要 求 , 到 达 喉 部 的 
声速 流 必 须 是 均匀 的 。 经 验证 明 , 如 果 稳 定 段 来 流 
是 均匀 的 ,只 要 有 一 条 光滑 连续 而 又 渐变 的 收缩 曲 
线 就 能 基本 满足 要 求 。 风 洞 试验 段 M 数 不 同 , 喉 部 
面积 也 不 同 , 因 而 收缩 比 是 随 W 数 而 变化 的 ,M 数 
越 高 ,收缩 比 越 大 。 亚 音速 段 外 形 采 用 移 轴 的 维 托 
辛 斯 基 公 式 , 并 由 喉 道 处 的 连续 性 条 件 来 确定 移 
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轴 量 。 
边界 层 修正 通常 采用 Sivells-Payne 方法 。 


1.2 超声 速 段 无 黏 位 流 型 面 计算 


1.2.1 喉 道 区 型 面 74 参数 计算 
假设 喉 道 区 为 一 元 等 精 流 ,经 验证 明 喉 道 区 型 
面 用 一 元 三 次 方程 能 满意 的 获得 要 求 的 气流 流动 状 
态 。 为 此 ,三 次 方程 如 下 。 
y= a0+ax+ Ax +ax (1) 
系数 oo aa as 由 边界 条 件 确定 ,边界 条 
件 ,为 


(2) 
加 由 式 (2) 代 入 式 (1) 经 整理 得 
Wo 下 ) > ? 
ul = 0 ? 
tang0， 
CQ2 全 » 
| 
tanO, 
Q3 二 一 3x2 (3) 
| 


“GG 以 式 (3) 代 入 式 (1) ,整理 后 得 
二 RE (二 _ a (二 i 
〇 式 (4) 为 所 求 的 喉 道 区 型 面 曲 线 ,其 中 0, 为 给 
定 的 经 验 值 。 x 由 下 述 关系 式 确定 ,为 
4 3 . ee © 
常规 高 超声 速 风 洞 ,气体 温度 在 1 000 K 左右 ， 
气体 仍然 可 以 作为 量 热 完全 气体 处 理 ,1 000 KK 以 
上 , 则 要 考虑 真实 气体 效应 。 
量 热 不 完全 的 真实 气体 喷 管 设计 方法 和 理想 气 
体 设计 方法 的 步骤 相同 ,但 原先 量 热 完全 的 气体 的 
速度 比 、 密 度 和 马赫 数 的 关系 不 再 能 适用 。 对 常规 
高 超声 速 风 洞 高 马赫 数 要 利用 NACA TR1135 的 修 
正 量 5 (同样 马赫 数 下 不 同 总 温 下 真实 气体 物理 量 
和 理想 气体 物理 量 的 比值 ) 进行 修正 。 
1.2.2 沿 径 向 流 波 前 4AB 气流 参数 计算 
在 波 前 4B 上 有 几何 分 割 关 系 ,为 


P 
Me = Ms - Wy(Ms -MD (6) 


式 中 : M4 < Ms < Ms 3;N 为 做 特征 线 网 格 时 沿 径 向 
流 波 前 48 的 分 割 总 数 , 从 妃 到 4,P=0,1,2,……,N。 
1.2.3 ABCD 转换 区 中 心 线 BC 上 的 导数 分 布 计算 

BC 上 的 导数 分 布 是 连续 变化 的 ,假定 速度 系 
数 成 三 次 方 分 布 


W = a0 +ax +aoX + aax (7) 
~ 和 一世 人 YY 
式 中 : % = 一 二 ;系数 a。、al 、w 、a3 由 边界 条 件 
C B 
决定 。 
在 B 点 
X=0,W=W,; 
| (ee) | 
dx 元 = 和 dx dx 计 守 dx X=Xp 
(8) 
在 C 点 
2 
2 
dx 5- dx js- 
(9) 


1.2.4 DCE 消 波 区 参数 计算 

利用 上 节 的 初始 条 件 ,做 特征 线 网 格 求解 型 面 
AD, 当 特征 线 网 格 做 到 DC, 则 以 DC 上 的 参数 作 初 
始 条 件 , 只 要 求 得 CE 上 的 初始 条 件 , 用 同样 的 方法 
作 特 征 线 网 格 ,就 可 以 解 出 型 面 DE, 于 是 整个 型 面 
AE 就 确定 了 。 

要 使 试验 段 获得 均匀 的 平行 气流 ,DCE 为 简单 
波 区 , CE 必须 为 一 条 直线 。 在 CE 上 的 M 数 均 为 
Ms 且 CE 上 的 参数 为 


Pp 
Xp=Xe 二 Je VME = 


Pp 
yp 一 N72 


(10) 


CQP | 
Ms 


Op =0 


1.3 轴 对 称 特 征 线 网 格 计 算法 


轴 对 称 喷 管 中 ,气流 特性 是 对 称 于 中 心 线 的 ,可 
以 只 研究 通过 中 心 轴 的 xy 平面 的 流 态 就 能 决定 整 
个 喷 管 的 流 态 。 

本 节 叙 述 用 轴 对 称 特征 线 网 格 计算 型 面 4D 和 
DE 的 方法 。 喷 管 的 对 称 轴 为 x 轴 ,气流 沿 x 轴 是 对 
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称 的 ,在 物理 平面 xy 上 设 已 知 点 己 (x,y ) 和 
P, (x%, ,2 ) ,PC ,yi) 的 左 特征 发 和 人 ,y ) 的 右 
特征 线 交 于 P (x; ,》3 ) ,如 图 2 所 示 。 


图 2 物理 平面 上 特征 线 与 参数 的 关系 


Fig.2 Relationship between feature lines and 


parameters on a physical plane 
1.4 ”壁面 坐标 计算 
FE 线 做 到 边界 上 的 点 4 后 ,为 了 要 确定 流 


点 发 出 的 左右 特征 线 。 
线性 插值 


和 一) 了 一 (11) 
到 X2 一 3 XR 一 %3 
得 
i (12) 
> < J 3 - 
6 由 流 线 和 特征 线 的 关系 
‘ A tang, 
二 MXR 一 XI (13) 
2 tan(C0 — a,) 
XR 一 %2 
由 上 式 解 出 
Re | es 
tangi (14) 


- tan0 - tan( 0b, — a,) 

{Cxitang -yi)-tangi[xotan(b — os) — 72]} 
(15) 

线性 关系 

XR—X 


“(ys 一 入 ) (16) 


NX3 一 %2 


yr = y+ 


A a a nw pe ns ee iNAax 
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7r。 AN 7A_ 出 口 十 
4 Iv 合作 日 口 十 
\IV OQ IFA) 


Mr = M+ Ms-M) (17) 
X3 一 2 

On 二 本) (18 ) 
Xa 一 2 


xx 由 式 (18) 给 出 。 式 中 :下 标 1 为 前 一 个 壁面 
点 ,初始 壁面 点 为 4 点 ;下 一 个 壁面 点 是 以 RR 为 初始 
条 件 进行 的 ,如 此 类 推 ,最 后 完成 4E 型 面 的 计算 。 


图 3 边界 上 流 线 与 特征 线 的 关系 


Fig.3 ”Relationship between streamlines and 


characteristic lines on the boundary 
2 超 音速 段位 流 型 面 边 界 层 修正 计算 


用 无 黏 流 特征 线 法 计算 得 到 的 喷 管 壁面 坐标 ， 
理论 上 应 该 就 是 喷 管 的 型 面 坐标 ， 但 实际 上 党 喷 管 
壁面 有 边界 层 增长 ， 傅 到 下 游 边界 层 愈 厚 。 边 界 层 
的 增长 ， 使 得 喷 管 气流 有 效 通道 截面 积 减 少 ， 因 而 
喷 管 出 口 不 能 得 到 设计 的 马赫 数 ， 并 且 也 不 能 得 到 
满意 的 均匀 流 场 。 为 防止 这 一 现象 ， 首 先 用 无 黏 特 
征 线 法 求 出 喷 管 无 黏 型 面 及 沿 壁面 马赫 数 分 布 ， 然 
后 把 它 作为 求 边界 层 位 移 厚 度 方 程 的 初始 值 ， 在 无 
黏 流 型 面 的 坐标 上 ， 加 上 边界 层 位 移 厚 度 ， 最 终 得 
到 喷 管 实际 型 面 坐标 。 采 用 这 种 方法 是 为 了 防止 因 
边界 层 的 存在 而 影响 试验 段 马 赫 数 的 分 布 的 均匀 性 
和 准确 性 。 高 超声 速 风 洞 喷 管 壁面 边界 层 一 般 ( 低 
密度 风 洞 除外 ) 为 满 流 边界 层 。 而 且 ， 在 风 洞 中 为 
了 防止 空气 组 分 液化 ， 需 要 把 空气 加 热 到 数 百度 至 
数 千 度 ， 因 此 ， 有 的 风 洞 喷 管 壁 面 还 需要 冷却 〈 常 
规 高 超声 速 风 洞 的 高 马赫 数 喷 管 ) 。 由 于 壁面 冷却 ， 
在 壁面 上 存在 热 交 换 ， 这 样 ， 喷 管 边界 层 计算 就 变 
得 复杂 。 

对 于 潮流 边界 层 的 近似 计算 方法 ,在 低速 时 ,有 
多 种 方法 ,而 用 于 总 温 很 高 的 高 超 音 速 的 近似 计算 
方法 一 般 有 4 种 , 即 Sivells-Payne 方法 , Reshtko- 
Tucker 方法 ,Persh-Lee 方法 和 Bartz 方法 。 边 界 层 
修正 计算 , 即 在 喷 管 的 每 一 个 横 截 面 上 给 以 一 个 质 
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量 补偿 量 , 它 对 超 音速 喷 管 也 同样 适用 ,不 同 之 处 在 
于 高 超 音 流动 需要 考虑 温度 效应 。 

对 有 热 交 换 和 压力 梯度 的 高 超声 速 边界 层 , 特 
别 是 高 烩 气体 边界 层 计算 方法 现在 还 不 太 成 熟 , 大 
多 数 是 把 低速 潮流 边界 层 计 算 方 法 移植 到 高 超声 速 
边界 层 , 用 得 比较 多 的 是 Sivells-Payne 方法 。 本 研 
究 采 用 Sivells-Payne 方法 进行 边界 层 修正 。 


2.1 动量 方程 的 建立 


Sivells-Payne 方法 先 用 修正 的 Stewartson 转换 ， 
简化 冯 . 卡门 轴 对 称 动量 积分 方程 ,在 解 方 程 时 ,用 
Eckert 参考 温度 法 及 气体 的 特性 ,用 修正 的 Crocco 
二 次 分 布 定律 给 出 边界 层 的 温度 分 布 ,这 样 即 把 方 
程 中 的 可 压缩 边界 层 形 状 系数 和 摩 阻 系数 变 为 不 可 


压缩 形式 ,又 考虑 了 非 绝 热 壁 温和 Pr 数 不 等 于 1 的 
影 租 。 


忆 ) 在 高 超 音 速 喷 管 中 ,边界 层 属于 灌流 边界 层 , 满 
流 壳 界 层 增 长 用 冯 “' 卡门 轴 对 称 动量 积分 方程 来 
2 求 边界 层 形状 因子 系数 内 

C J 采用 修正 的 Crocco 二 次 定律 ,可 给 出 考虑 了 绝 
热 睫 温和 Pr 关 1 影响 的 可 压缩 流 边 界 层 形状 因子 系 
数 昌 , ,边界 层 内 的 温度 分 布 为 

有 1 a 了 一 了 U 了 了 U 

‘ 心 了 ( 了 ) ( | 人 了 ) ( 


(19) 


式 审 : 7, 为 壁 温 ; 7,, 为 绝热 壁 温 。 


2.3 求 摩擦 系数 C， 


压缩 流 比 不 可 压缩 流 的 知识 不 够 完善 ,因此 ,使 
用 所 谓 参考 温度 的 概念 ,从 而 可 以 利用 不 可 压 流 的 
Cs 求 出 可 压 流 的 C1 ,所 谓 参 考 温度 法 ,就 是 即使 有 
压缩 性 存在 的 场合 ,在 壁面 附近 也 形成 不 可 压缩 流 ， 
所 以 把 适当 的 温度 作为 基准 来 评价 p 、y 。 那 么 ,不 
可 压缩 流 情况 下 的 关系 就 可 以 用 于 可 压缩 流 , 取 壁 
温 作为 参考 温度 时 ,可 以 知道 评价 压缩 性 大 小 程度 。 
利用 Eckert 给 出 的 参考 温度 法 , 即 参考 温度 了 为 


1 = 六 (7 + 7)+0.2207, - 7) (20) 


27 pi / ， 了 了 
全 二 人 - “= Cr 本 

pauses pr ps 7 
式 中 ,C= F(Re) 使 用 的 是 不 可 压 情 况 下 的 


Cr (21) 


八 \ (所 甘 口 工 
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Cs = 了 (Re,) 关系 式 。 


2.4 实际 计算 要 考虑 的 因素 


2.4.1 Re 数 的 范围 

Karmdn 一 Schooner 的 公式 ,只 适合 没有 压力 梯 
度 的 情形 ,因此 应 注意 到 式 (20) 的 结果 也 是 针对 没 
有 压力 梯度 的 。 当 logRe, = 1.5 或 logRe, = 2.3686 
时 , Cs 会 出 现 奇异 点 ,以 上 公式 适用 在 Re =10 ~10” 
范围 ,本 研究 的 计算 是 满足 要 求 的 。 
2.4.2 假想 原点 

边界 层 的 计算 通常 是 以 叭 道 为 初始 条 件 的 , 当 
x =0 时 , 求 得 C 是 发 散 的 ,为 了 消除 奇异 点 ,以 喉 部 
上 游 为 假想 原点 , 则 喉 部 处 x 短 0 ,假想 原点 的 x, 由 
下 式 求 得 。 


上 (22) 


和 三 


这 样 原来 的 x 坐标 从 喉 道 往 上 游 移动 x, 的 距 
离 , 移 轴 后 的 坐标 应 为 x,=x, +x 。 
2.4.3 边界 层 修 正 因子 

为 了 得 到 预先 给 定 的 出 口 直径 ,在 进行 边界 层 
修正 时 ,要 引入 边界 层 修正 因子 f ,修正 后 面 
坐标 为 


1 = hy, + 0. (x， ) (23 ) 


X, = fur 
计算 时 ,第 一 步 利 用 无 香 位 流 型 面 参数 ,计算 出 
喷 管 出 口 的 6" (x) ,然后 通过 式 (23) 计算 出 fi ,利用 
fh 进行 计算 ,得 到 修正 后 的 喷 管 型 面 参 数 YX ,再 
利用 XX 求 出 85 (x) 直到 满足 1 ys -71 
10mm 为 止 。 


3 ” 亚 音速 段 型 面 计算 


亚 音 速 段 型 面 计算 采用 维 托 辛 斯 基 公 式 , 高 内 
数 时 ,曲线 很 陡 ,必须 进行 适当 的 移 轴 , 才 能 较 好 满 
足 要 求 ,如 图 4 所 示 。 


图 4 亚 音速 段 型 面 坐标 


Fig.4 Subsonic section profile coordinates 
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为 了 保证 喉 道 处 气流 均匀 和 连续 ,在 该 处 应 满 
足 连续 条 件 , 即 收缩 段 喉 道 处 的 曲率 半径 应 等 于 边 
界 层 修正 后 的 超 音速 段 喉 道 处 的 曲率 半径 ,由 此 确 
定 移 轴 量 ,计算 按 以 下 公式 进行 。 


7 = , 3 (24) 
;| |1- 
Bal 


式 中 : ac =v37L ; 蕊 = 蕊 +D 
Y, =Y, +b;Y=Y -6b; 


1 _ 
I lla a 


n-l 
本 入 
> - 一 一 - (25) 
(© | 人 | 
LO n 


甚 中: Y 、 了. 和 了 分 别 为 收缩 段 的 人 口 , 喉 道 的 半 
息 和 工 条 上 的 截面 半径 ; 忆 、Y, 和 了 分别 为 移 轴 
后 收缩 段 的 入口 , 喉 道 的 半 高 度 和 工 轴 上 的 截面 半 
各 S 为 特征 长 度 ; b 为 对 X 的 移 轴 晤 ; n 表示 亚 音 
速 六 移 轴 后 的 进出 口 半 径 高 度 比 , 即 =Y/Y,。 


Co 
4 中 \ 轴 对 称 喷 管 设计 计算 程序 
> 


< 高 超声 速 风 洞 轴 对 称 喷 管 的 气动 设计 计算 程序 
使 用 MATLAB 软件 编制 ,利用 MATLAB 提供 的 丰富 
的 对 学 处 理工 具 函数 使 编程 相对 简单 。 程 序 主要 由 
5 部 分 组 成 :四 轴 对 称 喷 管 无 粘 流 型 面 计算 ;@ 边 界 
层 修正 计算 ;@ 收 缩 段 计 算 ; 四 结果 输出 整理 ;@ 相 
关子 函数 。 

依照 以 上 方法 , 表 1 给 出 了 M4.5 喷 管 型 面 设 
计 计 算 的 主要 结果 ,图 5 给 出 了 M4.5 的 型 面 曲线 。 
从 结果 可 以 看 出 , 当 过 了 C 点 后 ,马赫 数 就 达到 了 
4.5 ,表明 马赫 数 均 匀 区 较 长 ,可 以 满足 试验 要 求 。 

表 1 M4.5 无 黏 流 型 面 参数 主要 计算 结果 


Tab.1 The main calculation results of 


M4.5 non-viscous flow surface parameters 


Dot X/ mm Y/mm M W To0/(°) 
yh 0.000 61. 430 1.000 1.000 0.000 
4 926. 562 142.753 3.261 2.020 7.500 
B 1 596. 595 0.000 4.300 2.173 0.000 
CG 1 989.510 0.000 4.500 2.194 0.000 
E 3086.380 250.000 4.500 2.194 0.000 


~ ea eS AT 
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D 1 142.556 170.336 3.604 2.081 6. 823 


0 
-500 -200 100 400 700 10001300160019002200250028003100 
图 5 M4.5 喷 管 型 面 曲线 
Fig.5 M4.5 nozzle profile curve 


$5 M4.5 喷 管 流 场 分 布 特性 计算 分 析 


5.1 计算 方法 概述 


基于 我 们 目前 所 需要 的 结果 ,作为 验证 设计 程 
序 的 可 靠 性 来 讲 , 如 果 从 N-S 方程 出 发 ,通过 坐标 转 
换 、 方 程 离散 .边界 条 件 处 理 等 一 系列 复杂 数学 过 程 
到 最 后 编程 是 没有 必要 的 ,结合 实际 情况 ,利用 现 有 
的 CFD 软件 可 起 到 事半功倍 的 作用 。 

本 研究 对 M4.5 喷 管 流 场 分 布 特性 计算 采用 商 
用 Gambit 软件 来 生成 网 格 。GAMBIT 的 主要 优点 体 
现在 其 突出 的 非 结 构 化 的 网 格 生成 能 力 和 生成 过 程 
较 强 的 自动 化 能 力 ,包含 了 数 十 种 网 格 生成 方法 ,能 
够 针对 复杂 的 几何 外 形 生 成 三 维 四 面体 .六 面体 的 
非 结 构 化 网 格 及 混合 网 格 ,同时 Gambit 比较 注重 软 
件 的 图 形 用 户 界面 (GUI) ,交互 式 特点 直观 易学 。 
M4.5 喷 管 的 计算 网 格 如 图 6 所 示 , 网 格 内 部 间距 
10 mm。 当 网 格 生 成 完毕 ,要 定义 网 格 使 用 范围 ,并 
定义 好 边界 条 件 , 以 备 FLUNET 软件 计算 使 用 。 


图 6 M4.5 喷 管 流 场 计算 网 格 
Fig.6 M4.5 nozzle flow field calculation grid 


FLUENT 软件 提供 了 从 不 可 压 到 可 压 、 层 流 、 济 
流 等 范围 很 广 的 模拟 计算 能 力 。 软 件 中 输 运 现象 的 
数学 模型 与 所 模拟 的 几何 图 形 是 结合 在 一 起 的 。 满 
流 模型 是 FLUENT 中 很 重要 的 一 部 分 ,清流 会 影响 
到 其 它 的 物理 现象 如 浮力 和 可 压缩 性 。 满 流 模 型 具 
有 很 大 的 适用 范围 ,而 不 需要 对 特定 的 应 用 做 特殊 
调节 ,而 且 它 涵 括 了 其 它 物理 现象 的 影响 ,如 浮力 和 
可 压缩 性 。 通 过 使 用 扩展 壁面 函数 和 区 域 模型 ,对 
近 壁 面 精 度 问题 有 很 好 的 处 理 能 

本 研究 利用 FLUNET 软件 来 模拟 计算 前 面 设计 
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的 M4.5 喷 管 的 流 场 分 布 。 所 计算 喷 管 流 场 的 计算 
区 域 由 壁面 .压力 人 口 条 件 .压力 出 口 条 件 和 对 称 轴 
所 包围 的 区 域 组 成 。 因 本 研究 计算 喷 管 流 场 时 从 亚 
声速 到 超声 速 为 一 体 计 算 的 ,中 间 的 流速 状态 转换 
由 FLUNET 自动 模拟 , 故 计算 初 值 需要 多 次 设置 ,再 
加 上 喷 管 内 部 真实 的 边界 条 件 使 用 单个 的 计算 模型 
都 不 能 完全 真实 模拟 喷 管 内 部 的 流动 状态 , 故 计 算 过 
程 中 ,采用 不 可 压 ,C, 定 值 ,C, 变 值 ,理想 气体 逐步 计 
算出 的 。 下 面 给 出 了 采用 层 流 模型 计算 出 的 结果 。 

图 7 为 层 流 模型 计算 残 差 收 敛 过 程 图 ,图 中 开 
始 收敛 平稳 ,而 在 喷 管内 部 流动 由 汕 流 向 层 流转 换 
过 程 中 ,中 间 有 个 较 大 波动 ,但 很 快 又 平稳 收敛 。 从 
图 7 中 可 以 看 出 ,计算 收敛 过 程 较 快 ,说 明 网 格 质 量 
和 计算 所 选用 模型 较 好 。 
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图 7 计算 残 差 收敛 过 程 


g.7 The residual convergence process was calculated 


-一 对 于 高 超声 速 风 洞 来 说 速度 场 一 般 指 喷 管 内 和 
试验 段 的 马赫 数 分 布 规 律 , 本 文 主要 针对 速度 场 , 温 
度 场 , 方 向 场 进行 简要 分 析 ,目的 在 于 检验 本 文 计算 
的 喷 管 型 面 坐标 是 否 符合 设计 要 求 。 

M4.5 喷 管 的 马赫 数 分 布 如 图 8 所 示 ,从 图 中 可 
以 看 出 ,在 喷 管内 和 喷 管 出 口 处 的 马赫 数 分 布 趋势 
同 设计 分 布 趋 势 一 致 ,能 得 到 试验 所 需 的 流 场 参数 。 
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图 8 马赫 数 等 值 线 


Contours of Msch Number 


C hin aX iV 合 作 期 刊 | 第 40 卷 
Fig.8 Mach number contour 
在 喷 管 设计 方面 ,一 般 最 关心 的 是 喷 管内 部 和 
出 口 处 的 流 场 特性 情况 ,为 了 更 容易 比较 , 表 2 利用 
FLUENT 计算 结果 给 出 了 M4.5 喷 管 出 口 速度 场 指 
标 分 析 结 果 。 从 表 2 可 以 看 出 :在 最 大 马赫 数 偏差 
小 于 1.0% 的 情况 下 ,k-e 潮流 模型 计算 出 的 均匀 区 
直径 为 0. 37 m, 平 均 马 赫 数 为 4. 479, 均 方差 为 
0. 005 3 ;而 层 流 模 型 计算 出 的 均匀 区 直径 为 0.48 m， 
平均 马赫 数 为 4.499, 均 方差 为 0.0090; 这 满足 参考 
文献 15] 中 规定 的 指标 要 求 ,同时 也 符合 参考 文献 
[6] 中 规定 的 指标 要 求 。 
表 2 M4.5 喷 管 出 口 速度 场 分 析 结 果 
Tab.2 Analysis results of exit velocity field of M4.5 nozzle 


计算 模型 清流 i-e 层 流 1 层 流 2 
D 0.37 0.48 0.42 
MavE 4.4794 4.4987 4.5008 
有 4.4853 4.509 2 4.509 2 
JiN 4.443 6 4.4753 4.4893 
Ory 0.0053 0.0090 0.0065 
DM,_yax 0.0357 0.023 4 0.0115 


而 从 高 超声 速 风 洞 现 有 的 喷 管 流 场 校 测 结果 来 
推断 ,实际 加 工 出 的 M4.5 喷 管 的 流 场 均匀 区 直径 
应 在 两 者 之 间 ,估计 0.42 m 是 比较 合理 的 , 层 流 模 
型 计算 出 的 结果 取 均 匀 区 直径 为 0.42 m 时 ,平均 马 
赫 数 为 4. 500 , 均 方 差 为 0.006 5; 这 满足 参考 文献 
[5] 中 规定 的 先进 指标 要 求 , 同 时 也 符合 参考 文献 
[6] 中 规定 的 指标 要 求 。 可 见 如 果 喷 管 加 工 按 要 求 
进行 ,本 套 M4.5 喷 管 型 面 数据 是 符合 要 求 的 。 


6 结 论 


通过 对 高 超声 速 风 洞 马赫 数 4.5 喷 管 的 内 型 面 
计算 ,并 对 流 场 进 行 计算 分 析 , 可 以 得 出 以 下 结论 。 

1) 本 研究 采用 的 轴 对 称 喷 管 气动 设计 计算 方法 
得 出 的 喷 管 马赫 数 与 所 需 马赫 数 误差 在 1.0% 以 
内 ,满足 国 军 标 对 喷 管 马赫 数 的 设计 要 求 。 

2) 本 研究 在 数值 验证 计算 中 ,采用 层 流 模型 计 
算 残 差 ,其 中 连续 方程 .x 和 y 方向 速度 ,能 量 方 程 的 
残 差 都 衰减 到 10“ ~10“, 收 敛 很 好 ,说 明 采 用 的 计 
算 模型 合理 。 

3) 本 研究 通过 马赫 数 4.5 喷 管 的 气动 设计 和 数 
值 验证 ,得 出 了 满足 国 军 标 要 求 的 喷 管内 型 面 数据 ， 


投稿 网 站 :http://ejam. xjtu. edu. en ” 微 信 公众 号 :应 用 力学 学 报 


I 
第 2 其 |， 等 ， 高 超声 速 风 洞 马赫 数 4 5 喷 管 气动 设计 和 数值 的 尼 naXiv 合 作 期 和 


黄 风 ， 等 


可 用 于 M4.5 喷 管 的 设计 加 工 。 
参考 文献 : 


[1] 


MD 


纺 


chinaXiv 


202307.00056v 


黄 确 , 杨 永 能 , 刘 奇 ,等 . 高 超声 速 风 洞 马赫 数 4.5 喷 
应 用 [可 .西北 工业 大 学 学 报 ,2021 ,39(5) :1064-1069. 
HUANG Ju, YANG Yongneng,LIU Qi,et al. Developing and apply- 


管 研 制 与 


ing Mach 4. 5 nozzle in hypersonic wind tunnel [ J ]. Journal of 
Northwestern Polytechnical University ,2021 ,39(5) :1064-1069 (in 
Chinese ) . 

中 国人 民 和 解放军 总 装备 部 . 高 超声 速 风 洞 气动 力 试验 方法 : 
GJB 4399 一 2002[ Sj. 北京: 中国 人民 解放 军 总 装备 部 ,2002. 
周 晓 刚 , 白 本 奇 . FL-31 风 洞 岗位 操作 指南 [2Z]. CARDC- 
2 ,2007. 

杨 永 能 , 黄 刚 . FL-31 风 洞 M4. 5 喷 管 
空气 动力 研究 与 发 展 中 心 ,2019. 
张 饥 钊 . PL31 风 洞 号 管 气动 设计 计算 方法 [R]. 中 国 空气 动 
力 研究 与 发 展 中 心 高 速 空气 动力 研究 所 ,1976. 

HAGEMANN C,IMMICH H. Critical assessment of the linear plug 


型 面 设计 计算 [R]. 中 国 


加 


2020 


[10] 


西部 科技 期刊 联盟 


BS 


309 


nozzle concept[ C ]//37th Joint Propulsion Conference and Exhib- 
it. Reston, VA , USA :AIAA ,2001 :AIAA 2001-3683. 
ROWBOTHAM M. XRS-2200 linear aerospike engine - use of Pro/ 
ENGINEER for determining mass properties[ C |]//35th Joint Pro- 
pulsion Conference and Exhibit. Reston, VA, USA: AIAA, 1999 : 
AIAA 1999-2334. 

BOOTH T E, VILJA J O,CAP D P,et al. The design of linear aero- 
spike thrust cells[ C ]//29th Joint Propulsion Conference and Ex- 
hibit. Reston, VA , USA : AIAA ,1993:. AIAA 1993-2562. 
MEENAKSH R, HOFFMAN J D, MURTHY S N B. Design and 
performance computations in complex 3-D nozzles[ C |]//37th Aero- 
space Sciences Meeting and Exhibit. Reston, VA, USA.: AIAA, 
1999 :AIAA 1999-882. 
VAN LEER B. Towards the ultimate conservative difference 
scheme. V. A second-order sequel to Godunov’s method[ J ]. Journal 
of computational physics ,1979 ,32(1 ) :101-136. 


(编辑 史 淑 英 ) 


应 用 力学 学 报 


被 评 为 2022 年 度 “ 西 牛 计划 ”之 


精品 中 文科 技 期 刊 


投稿 网 站 :http://cjam. xjtu. edu. cn 


微 信 公 众 号 :应 用 力学 学 报 


